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The recent progress in understanding the QCD phase diagram and the physics of the 
QCD critical point is reviewed. 

§1. Preamble 

Quantum Chromodynamics is one of the most remarkable theories of Nature. 
Its mathematical foundations are concise, yet the phenomenology which the theory 
describes is broad and diverse. QCD phenomenology at finite temperature and 
baryon number density is one of the least explored regimes of the theory. There 
are several experimental windows into this regime. One is the physics associated 
with the interior of neutron stars. Another, which is the subject of the ongoing and 
planned experimental programs, is the physics of heavy ion collisions. 

This report is a review of recent developments in our understanding, mostly 
theoretical, but also experimental, of the phase diagram of QCD. There exist a 
number of excellent recent reviews 1 ^ 6 - 1 which discuss many of the questions addressed 
here as well as the related material not covered in this report. The report focuses 
on the physics of the critical point of QCD and its search. There are many open 
questions in this field, and some of the theoretical as well as experimental results and 
expectations discussed here might not hold under further scrutiny. Nevertheless it 
is hoped that this report will provide a useful contemporary guide to both theorists 
and experimentalists entering the field as well as a stimulating reading to the field's 
experts. 

§2. What is the QCD critical point? 



2.1. The phase diagram 

Fig. 1 shows a sketch of the QCD phase diagram as it is perceived by a modern 
theorist. By a phase diagram we shall mean the information about the location of 
the phase boundaries (phase transitions) as well as the physics of the phases that 
these transitions delineate. The phase transitions are the thermodynamic singular- 
ities of the system. The system under consideration is a region (in theory, infinite) 
occupied by strongly interacting matter, described by QCD, in thermal and chemical 
equilibrium, characterized by the given values of temperature T and baryo-chemical 
potential In practice, it can be a region in the interior of a neutron star, or 
inside the hot and dense fireball created by a heavy ion collision. 
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Fig. 1. QCD phase diagram 



On the phase diagram, the regime of small T and large ^b is of relevance to neu- 
tron star physics. Because of low temperature, a very rich spectrum of possibilities 
of ordering can be envisaged. The line separating the Color-Flavor-Locked (CFL) 
phase, predicted in Ref. 7), from the higher temperature disordered phase (quark- 
gluon plasma, or QGP) is the most simplified representation of the possible phase 
structure in this region. This regime is also of particular theoretical interest because 
analytical controllable calculations are possible, due to asymptotic freedom of QCD. 
The reader is referred to the reviews 1 ) -5 ) which cover the recent developments in the 
study of this domain of the phase diagram. 

The region of the phase diagram more readily probed by the heavy ion collision 
experiments is that of rather large T ~ 100 MeV, commensurate with the inherent 
dynamical scale in QCD, and small to medium chemical potential [lb ~ — 600 
MeV. Theorists expect that this region has an interesting feature - the end point 
of the first order phase transition line, the critical point marked E on Fig. 1. The 
physics of this point is the focus of the review. 

2.2. Why should there be a critical point? 

The argument (which is not a proof) that the point E must exist is short, and 
is based on a small number of reasonable assumptions. The two basic facts that it 
relies on are as follows: 

(1) The temperature driven transition at zero \xb is not a thermodynamic sin- 
gularity. Rather, it is a rapid, but smooth, crossover from the regime describable 
as a gas of hadrons, to the one dominated by internal degrees of freedom of QCD - 
quarks and gluons. This is the result of finite T lattice calculations. 8 ) 

(2) The hb driven transition at zero T is a first order phase transition. This 
conclusion is less robust, since the first principle lattice calculations are not con- 
trollable in this regime (naive Euclidean formulation of the theory suffers from the 
notorious sign problem at any finite hb)- Nevertheless a number of different model 
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approaches 9 ) -16 ) (see Section 3.3) indicate that the transition in this region is strongly- 
first order. 

(3) The last step of the argument is a logical product of (1) and (2). Since the 
first order line originating at zero T cannot end at the vertical axis \ib = (by virtue 
of (1)), the line must end somewhere in the midst of the phase diagram. 

The end point of a first order line is a critical point of the second order. This 
is by far the most common critical phenomenon in condensed matter physics. Most 
liquids possess such a singularity, including water. The line which we know as the 
water boiling transition ends at pressure p = 218 atm and T = 374°C. Along this 
line the two coexisting phases (water and vapor) become less and less distinct as one 
approaches the end point (the density of water decreases and of vapor increases), 
resulting in a single phase at this point and beyond. 

In QCD the two coexisting phases are hadron gas (lower T), and quark-gluon 
plasma (higher T). What distinguishes the two phases? As in the case of water and 
vapor, the distinction is only quantitative, and more obviously so as we approach the 
critical point. Rigorously, there is no good order parameter which could distinguish 
the two phases qualitatively. The chiral condensate, (tpip), which comes closest to 
being an order parameter, is non-zero in both phases because of the finite bare 
quark mass. Deconfinement, although a useful concept to discuss the transition 
from hadron to quark-gluon plasma, strictly speaking, does not provide a good order 
parameter. Even in vacuum (T = 0) the confining potential cannot rise infinitely - a 
quark- antiquark pair inserted into the color flux tube breaks it. The energy required 
to separate two test color charges from each other is finite if there are light quarks. 

2.3. Critical or tricritical? 

There is an idealization of QCD where the distinction between the hadron gas 
and quark gluon plasma is sharp. It describes the world with massless quarks m q = 0. 
In this limit of QCD with 2 massless quarks (up and down) the chiral symmetry 
SU(2)y xSU(2)a is exact. Although interactions respect this symmetry, it is sponta- 
neously broken in the QCD vacuum to SU(2)y , and the Goldstone theorem demands 
3 massless bosons - the pions. This breaking is a result of nonperturbative dynamics 
in QCD (instantons provide a natural mechanism). 

The breaking of the global symmetry, such as the chiral symmetry, can be 
thought of as establishment of the long-range order in the vacuum. It is the or- 
der which dictates the preferred SU(2)^ direction for all points in space, and over 
which the pions are quantized "ripples" . At sufficiently high T, the order is melted 
as in any other such system (compare, e.g., to the disordering of the ferromagnet 
at Curie temperature). The chiral symmetry is restored. The two phases must be 
separated by a thermodynamic singularity - a phase transition. 

This argument can be made more rigorous by considering the order parameter, 
(V>V')> the expectation value, or the condensate, of a field transforming non-trivially 
under the broken symmetry. Consider a fixed value of hb, e.g., (1b = 0. The 
chiral condensate (ipip) as a function of T is identically zero (by symmetry) for all 
temperatures above some value T c and is nonzero function of T below (symmetry 
breaking). Such a function cannot be analytic. A singularity must occur at T c . 
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In QCD, lattice calculations show that this singularity is a second order phase 
transition if T c is approached at \i& = 0. 8 ^' 17 ^ 23 ^ At other values of hb the critical 
temperature T c is different, but the line of transitions T c (/ib) cannot terminate, since 
any path from the vacuum T = \xb = to the high T phase must cross a singularity. 
Somewhere in the midst of the phase diagram the order of the transition should 
change to first order (according to point (2) of Section 2.2).*) The point where this 
happens is the tricritical point. The resulting phase diagram is illustrated in Fig. 2. 




Fig. 2. Phase diagram of QCD with two massless quarks. The chiral symmetry order parameter 
qualitatively distinguishes two phases: {i/>V) 7^ in the broken phase and (4>il>) = in the 
symmetric phase. 



Once the quark mass m q is turned 
back on, the distinction between the 
symmetric and broken phases is blurred, 
and the second order phase transition 
is replaced by a smooth crossover. The 
situation is analogous to the ferromag- 
net — an arbitrary small magnetic field 
(the analog of m q ) smooths away the 
Curie singularity (Fig. 3). The first or- 
der phase transition, on the other hand, 
is associated with a finite discontinuity 
of the order parameter and cannot be 
removed by an arbitrarily small pertur- 
bation m q 7^ 0. Thus we arrive back at 
the diagram in Fig. 1. 

It is also useful to take a look at Fig. 



M 




T c T 



Fig. 3. The order parameter vs temperature in 
a Curie ferromagnet with zero and non-zero 
applied magnetic field. In QCD, the chiral 
order parameter {'tpip} behaves similarly as 
a function of T at m q = and m q 7^ 0. 

4, where the 2-dimensional T[±b phase 



We also assume that there is only one transition between the broken and symmetric phases. 
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Fig. 4. A three-dimensional view (T, /j, B , m q ) of the QCD phase diagram near the tricritical point. 

diagram is extended to 3-dimensions by adding the quark mass m q as the third axis. 
One can see that the second order transition line at m q = does not extend into 
m q 7^ 0. This line can be seen as a boundary of the coexistence surface of the two 
spontaneously broken phases with (ipip) of opposite signs. A first order line ending 
at a critical point, on the other hand, exists for all nonzero (small) m q , thus making 
up a surface which looks like two wings in Fig. 4. The tricritical point can be seen 
as the end of a first order line where 3 phases coexist (line of triple points). 

Another useful sketch is made in Fig 5. It shows, in a schematic way, the shape 
of the effective potential in various regions around the tricritical point. One can 
see that the three minima, which are equally deep on the triple line fuse into one 
minimum at the tricritical point. 




(triple line 

Fig. 5. Illustration of the shape of the effective potential for the chiral order parameter near the 
tricritical point in the m q — plane. Two additional (spinodal) lines, not present in Fig. 4 
indicate the boundary of the existence of metastable minima. 
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2.4. Critical behavior: static and dynamic universality class 

Determining properties of QCD (equation of state, correlation functions, etc.) 
near the critical point is difficult, for the same reason as it is difficult to find the 
location of the critical point (see next Section). However, as it is the case for any 
critical point, singular properties, such as critical exponents, can be determined using 
universality arguments. 

According to the scaling postulate, central to the theory of critical phenom- 
ena, 24 ) all singular contributions to the thermodynamic quantities are powers of the 
correlation length £, which diverges at the critical point. These powers, or critical ex- 
ponents, are universal, in the sense that they depend only on the degrees of freedom 
in the theory and their symmetry, but not on the other details of the interactions. 
Very different physical systems may belong to the same universality class, as far as 
their critical behavior is concerned. 

One should distinguish static and dynamic universality classifications. 25 ' 1 From 
the point of view of static critical phenomena, the QCD critical point falls into the 
universality class of the Ising model. This is a consequence of the fact that at m q 7^ 
no symmetry remains which would require the order parameter to have more than 
just one component. The field theory which describes the static critical behavior, 
the one-component </> 4 theory in 3 dimensions, has the critical exponents of the Ising 
model. *' 

What is the nature of this order parameter? It can be taken as the value of the 
chiral order parameter tpip (often called a), since it is distinct in two phases coex- 
isting across the first order phase transition terminating in the critical point. As a 
result, the static (equal-time) correlation function (vpijj(x)if)ip(y)) develops divergent 
correlation length: 

1 



|as|»£; 



(2-1) 



where (^(x)ifjif)(0)} c = (^(x^tp^)) — (4>?p} 2 . The correlation length diverges, 
£ — > 00, as we approach the critical point. For the Ising universality class r] ~ 0.04. 

Another interesting quantity, both from theoretical and experimental points of 
view, is the baryon number density ns{x). Because symmetry (or, rather, the ab- 
sence of such) allows mixing of ns{x) with vpvp(x), the divergence of the baryon 
number susceptibility is related to the divergence of the correlation length £: 



dn B 
dfi B 



J d 3 x(n B (x)n B (0)) c ~ J d 3 x(i;i;(x)^(0)} c ~ f^. (2-2) 



The baryon number density also jumps across the first order phase transition. 
One can equally well use ub as the degree of freedom in the effective theory near 

As another example, consider any of the critical points on the 2nd order line at m q = on 
Fig. 2. Because of the 0(4)~ SU(2)y xSU(2)a symmetry, which is restored at this critical point, the 
order parameter must carry 4 components — sigma and 3 pions (a, n). The resulting field theory 
describes the universality class of the 0(4) ferromagnet. 26 ' -28 ' 
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the critical point, or any linear combination of V>V and iib (or any other field which 
can mix with ipi/j) which is discontinuous across the first order phase transition. 
Regardless of the choice, there is only one order parameter, as far as the static 
critical behavior is concerned. 

The situation is similar but a little more complicated if one considers dynamic 
critical behavior, e.g., the singularities of kinetic coefficients, etc. The scaling pos- 
tulate is similar in this case, but the universality classes are now determined by the 
degrees of freedom which define the effective hydrodynamic theory near the critical 
point. 25 ) In this case the fundamental difference between ipi/j and hb fields is that 
the latter is a conserved density. The hydrodynamic equations for ub are diffusive, 
while the dynamics of ipip is relaxational. Because the two modes mix, there is, again, 
only one independent hydrodynamic variable, and it is diffusive. 29 )' 30 ) This mode 
involves fluctuations of both ipip and ub in a fixed proportion. The fluctuations of 
•t/iip alone relax on a finite time scale even at the critical point.*) 

The complete hydrodynamic theory near the critical point must also involve 
the energy and momentum densities. Once the hydrodynamic equations are written 
down, and the mixing between vp-ip, ub and the energy density is taken into account, 
one finds the theory equivalent to the one describing the liquid-gas phase transition, 
model H in the classification of Ref. 25). One consequence of this theory, interesting 
from phenomenological point of view, is the vanishing of the baryon number diffusion 
rate at the critical point: D ~ (,~ XD , with exponent xd ~ l- 30 -* 

§3. Where is the critical point? Theory 

Theoretically, finding the coordinates (T, of the critical point is a well- 
defined task. We need to calculate the partition function of QCD and find the 
singularity corresponding to the end of the first order transition line. The Lagrangian 
of QCD is known, and the partition function is given by a path integral of the expo- 
nent of the QCD action, after Wick rotation to the Euclidean space (with imaginary 
time compactified on a torus of circumference 1/T). 

Of course, calculating such an infinitely dimensional integral analytically is be- 
yond our present abilities (perturbation theory is not an option here, in the relevant 
region of T and ub)- We are thus left with numerical methods, i.e., lattice Monte 
Carlo simulation. At zero ub this method allows us to determine the equation of 
state of QCD as a function of T and reach the conclusion (1) in Section 2.2. However, 
at finite \ib the Nature guards its secrets better. 

3.1. Importance sampling and the sign problem 

The notorious sign problem has been known to lattice experts since early days 
of this field. Calculating the partition function using Monte Carlo method hinges on 
the fact that the exponent of the Euclidean action Se is a positive definite function of 
its variables (values of the fields on the lattice). This allows one to limit calculation 



*' A related observation, that the sigma pole mass does not vanish at the critical point in the 
large- iV NJL model, was made in Ref. 14) and confirmed in Ref. 29). 
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to a relatively small set of field configurations randomly picked with probability 
proportional to the value of exp(— Se)- The number of such configurations needed 
to achieve reasonable accuracy is vastly smaller than the total number of possible 
configurations. The latter is exponentially large in the size V of the system, or, 
the number of the degrees of freedom: exp(const • V). The method, also known as 
importance sampling, utilizes the fact that the vast majority of these configurations 
contribute a tiny fraction because of the exponential suppression by exp(— Se)- Only 
configurations with sizable exp(— Se) are important. 

At fi 7^ the action Se is complex. What configurations are important then? 
A number of ways to circumvent the problem have been tried. For example, using 
the modulus of exp(— Se) as a measure of importance, or the value of exp(— Se) at 
zero fj,B- Unfortunately, none work, at present. 

3.2. The overlap problem 

For the latter choice, exp(— Se)\(i b =o, the problem can be understood physi- 
cally and is known as the overlap problem. The important configurations at \xb = 
are different than those of /xb ^ 0. How bad is this quantitatively? At finite vol- 
ume, even at (is = 0, the configurations important for \ib ^ pop up, but with 
a small probability. This probability is exponentially small as volume V — ► oo: 
exp(— const • V). When we calculate the partition function using this method, we 
correct for this by multiplying the contribution of these rare configurations by the 
factor exp(+5£;| Ms= o — Se)- The procedure is termed reweighting.*^ The reweight- 
ing factor is exponentially large as V — ► oo — both the magnitude and the complex 
phase are exp (const • V). Fluctuations, or statistical noise, in the exponentially tiny 
number of the rare important configurations completely washes out the significance 
of the result. 

In layman's terms, imagine that we want to study ice, but can only run exper- 
iments at normal room temperature and pressure. Using the reweighting method 
is analogous to trying to glimpse the information by waiting for rare configurations 
when all the water molecules accidentally gather in one corner of the lab, forming a 
chunk of ice. The amount of time that this experiment would require is exponentially 
large as V — > oo. 

3.3. Theoretical predictions 

The first lattice prediction for the location of the critical point has been reported 
in Ref. 32). The assumption is that, although the problem becomes exponentially 
difficult as V — ► oo, in practice, once can get a sensible approximation at finite V. In 
addition, simulations at finite T might suffer lesser overlap problem because of large 
thermal fluctuations. 33 ) One can hope that if the critical point is at a small value 
of fiB, the volume V may not need to be too large to achieve a reasonable accuracy. 
In particular, numerical estimates show 34 ) that the maximal value of he which one 
can reach within the same accuracy shrinks only as a power of l/V. 

However, it is not possible to determine this accuracy, since the exact result is 

The reweighting method in application to finite /j,b calculations is known as the "Glasgow 
method" (reviewed in Ref. 31)). 
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unknown. Normally, one would estimate the error by going to increasingly large 
volumes V, but, as discussed above, the method becomes prohibitive too quickly 
(exponentially) in this limit. Ultimately, the result of Ref. 32) might turn out to 
be a good approximation to the exact answer, but we can only tell once we have 
an independent result to compare it to. A qualitatively new approach is needed to 
overcome the QCD sign problem. *) 
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Fig. 6. Theoretical (models and lattice) predictions for the location of the critical point. The 
labels correspond to Table I. The two dashed lines indicate the magnitude of the slope d 2 T/dfi 2 
obtained by lattice Taylor expansion. 38 ' The upper curve agrees with Ref. 39). The lower curve 
corresponds to smaller quark mass. Errors/uncertainties are not shown. 



In the absence of a controllable (i.e., systematically improvable) method, one 
turns to model calculations. Many such calculations have been done. 9 ) -16 ) Figure 6 
and Table I summarize the results. One can see that the predictions vary wildly. An 
interesting point to keep in mind is that each of these models is tuned to reproduce 
vacuum, T = (j,b = 0, phenomenology. Nevertheless, extrapolation to nonzero /xb is 
not constrained significantly by this. In a loose sense, the existing lattice methods 
can be also viewed as extrapolations from \xb = 0, but finite T. 

Two new lattice approaches are being developed. Each of them has the capacity 
to determine the location of the critical point. One approach is based on simula- 
tions at finite imaginary values of fiB 39 ' and the other on Taylor expansions around 
fj,B = 0. 38 ) The curvature of the phase transition line found using these methods is 
indicated by the upper parabola in Fig. 6. Recent result 38 ) (lower parabola in Fig. 6) 
seems to indicate large sensitivity of this curvature to the quark mass. This may or 
may not be related to the strong sensitivity of the position of the critical point to 
the value of the strange quark mass observed in Ref. 39). Qualitatively, one should 

*^ In theories similar, or approximating, the finite density QCD, the sign and/or overlap prob- 
lems have been tackled recently, using various new methods see, e.g., Refs. 35)-37). 
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expect that the critical point moves toward smaller /ig as the strange quark mass 
m s is decreased, since for sufficiently small m s the chiral transition must turn first 
order according to renormalization group arguments. 26 ^ 



Source 


(T,(jlb), MeV 


Comments 


Label 


MIT Bag/QGP 


none 


only 1st order, no chiral symmetry 




Asakawa,Yazaki '89 


(40, 1050) 


NJL, CASE I 


NJL/I 




(55, 1440) 


NJL, CASE II 


NJL/II 


Barducci, et al '89-94 


(75, 273)tcp 


composite operator 


CO 


Berges, Rajagopal '98 


(101, 633)tcp 


instanton NJL 


NJL/inst 


Halasz, et al '98 


(120, 700)tcp 


random matrix 


RM 


Scavenius, et al '01 


(93,645) 


linear cr-model 


LSM 




(46,996) 


NJL 


NJL 


Fodor, Katz '01 


(160, 725) 


lattice reweighting 




Hatta, Ikeda, '02 


(95, 837) 


effective potential (CJT) 


CJT 


Antoniou, Kapoyannis '02 


(171, 385) 


hadronic bootstrap 


HB 


Ejiri, et al '03 


(?,420) 


lattice Taylor expansion 





Table I. Theoretical predictions of the location of the critical point. The predictions for tricritical 
point are indicated as 'TCP'. The last column gives the corresponding label on Fig. 6. 



§4. Scanning QCD phase diagram 

Even though the exact location of the critical point is not known to us yet, the 
available theoretical estimates strongly indicate that the point is within the region 
of the phase diagram probed by the heavy-ion collision experiments. This raises the 
possibility to discover this point in such experiments. 40 ' 
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Fig. 7. Example trajectories traced by a fireball created in a heavy ion collision on the phase 
diagram of Fig.l. Increasing the collision energy one moves the freezeout point (empty circle) 
to smaller /i_g, approaching and then receding away from the critical point. 

The idea is illustrated in Fig. 7. It is known empirically that with increasing 
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collision energy, yfs, the resulting fireballs tend to freezeout at decreasing values of 
the chemical potential. This is easy to understand, since the amount of generated 
entropy (heat) grows with yfs while the net baryon number is limited by that number 
in the initial nuclei. 

The trajectories on Fig. 7 terminating in the freezeout points indicate (theoreti- 
cally perceived) time history of a small, but thermodynamically macroscopic, volume 
of the expanding fireball, at various initial collision energies. In the approximation 
of ideal hydrodynamics these trajectories follow lines of constant baryon per en- 
tropy ratio (baryon asymmetry), due to conservation of the baryon number and the 
entropy. The characteristic discontinuity of the trajectory at the first order phase 
transition is a result of the discontinuity of the baryon asymmetry across this tran- 
sition. Because the shift is toward the end point E, this leads to the phenomenon 
of focusing: 40 ^ the freezeout points tend to cluster near the critical point for a wide 
range of initial trajectory points. Therefore, wider range of leads to freezeout in 
the critical region, making the task of finding the point somewhat easier. 
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Fig. 8. The freezeout points for different heavy-ion collision experiments. On the right, the lattice 
results 32 - 1 ' 38 "'' 39 '' are superimposed. Figures are reproduced from Ref. 41). 



The information about the location of the freezeout point for given experimental 
conditions is obtained by measuring the ratios of particle yields (e.g., baryons or 
antibaryons to pions), and fitting to a statistical model with T and /jg as parameters. 
Such fits are amazingly good, 41 ) and the resulting points for different experiments 
are shown in Fig. 8. 

For comparison, the location of the lattice reweighting calculation result 32 ) is 
superimposed on the experimental freezeout curve in Figure 8(right). It may appear 
that the critical point is somewhat away (higher T) from the freezeout curve. How- 
ever, as emphasized already in Section 3.3, the systematic error of the lattice result 
is not known, since the volume V — > oo limit is unattainable using the reweighting 
method. Even if one takes the lattice result at face value, one still has to take into 
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account the fact that the position of the critical point must shift to smaller values 
of fj,B once the quark masses (notably, the strange quark mass) are reduced toward 
their true values, from those used in Ref. 32). More recent lattice results 38 ^' 39 ) sup- 
port this expectation. The results of the model calculations (see Figure 6 and Table 
I) are closer to the freezeout curve than the band in Fig. 8. 

Additional effect, which plays a significant role, is the critical slowing down near 
the point i£. 30 )' 42 )' 43 ) This phenomenon is the major limiting factor (size limitation 
is less stringent 42 -*) for the maximal correlation length that can be achieved realis- 
tically in a heavy ion collision experiment. Although, as a result, the sharpness of 
the signatures of the critical point is reduced, another consequence is the shift of 
the position (due to delay) of the maximum of the correlation length toward lower 
temperatures . 43 ) 

On the experimental side, one has 
some control over the freezeout temper- 
ature by adjusting the size of the ions. 
Smaller systems freeze out somewhat 
earlier (higher T). Also, due to the sin- 
gularity in the specific heat, the freeze- 
out occurs at higher T for trajectories 
passing near the critical point. 40 ) 

Another interesting observation, 
potentially important for the search of 
the critical point is that of the shape of 
the critical region. 16 ) One can see in 
Figure 9 a model calculation of a di- 
vergent susceptibility near the critical 
point, which shows that the critical re- 
gion is stretched in the direction of the crossover line. This shape is easy to under- 
stand remembering that, in the (T, //B,m 9 ) space (see Fig. 4), the critical point is 
connected to the tricritical point (the black dot in Fig. 9) by a whole line of critical 
points. 

§5. Signatures: event-by-event fluctuations 
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Fig. 9. The model calculation 16 ' of the shape 
of the critical region. Note that fi = (ib/3. 



One of the actively pursued signatures of the critical point is the non-monotonous 
dependence on y/s (and thus, on (jlb) of the event-by-event fluctuation observ- 
ables. 40 )' 42 ) The idea can be understood qualitatively by noting that: (1) the suscep- 
tibilities diverge at the critical point, and (2) the magnitude of the fluctuations are 
proportional to the corresponding susceptibilities. For example, for the fluctuations 
of energy or charge, the well-known relations are 



(5-1) 



Ideally, one could determine susceptibilities on the left-hand sides by measuring 
the fluctuations on the right-hand side. 44 ) However, practically, the measurement of 
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the corresponding fluctuations, AE or AQ, is not feasible because not all the particles 
end up in the detector. 42 )' 45 ) A more differential measure of the fluctuations needs 
to be computed in theory and compared to experiment. 

5.1. Two-particle correlator 

A number of such measures can be obtained starting from a two particle corre- 
lator 

(An«AnP k ) = (n^ k )-(n«)(4) (5-2) 

where An® = n." — (n") is the event-by-event fluctuation of the number of particles of 
the type a in the momentum bin centered around p. Experts familiar with Hanbury- 
Brown-Twiss (HBT) interferometry 46 ) may recognize in (5-2) the HBT correlation 
function. 

The two-particle correlator (5-2) can be directly measured. However, for such 
a function of many variables, it might be difficult to represent the result of this 
measurement. A useful representation, for example, is obtained by limiting (pro- 
jecting) the variables to transverse components of p and k. The resulting plot of a 
function of two arguments, px and for, is often referred to as a 'Trainor plot' (see, 
e.g., Ref. 47)). Interesting information can be also obtained by projecting onto the 
rapidities of p and k. If in addition, one weights each particle with its charge, the 
resulting correlator, as a function of the rapidity difference y p — yk, is essentially the 
balance function introduced in. 48 ) 

There also exist many cumulative measures, proposed by theorists and/or used 
by experimentalists, 42 )' 45 )' 49 ) -55 ) that can be expressed in terms of correlator (5-2). 
As an example, the fluctuation of electric charge is given by summing over momenta 
p and k of all particles in the experimental acceptance window and weighting each 
particle with its charge q a : 

AQ = ^2q a An«; thus ((AQ) 2 ) = £ £ q^(An a p Ani). (5-3) 

P,a p,a k,/3 

The same applies to the fluctuations of the baryon charge, with q a substituted by 
the baryon charge of the particles. Similar equation (see Eq.(5-5)) also applies to 
the fluctuations of the mean transverse momentum px, in which case q a should 
be replaced with pr — pr - the deviation of the momentum px from the all-event 
(inclusive) mean pf. 

The correlator (5-2) can, in principle, be calculated, under assumption of thermal 
equilibrium, once the relevant interactions are known. In the case of the critical 
point, we need to concern ourselves with the interactions which can lead to singular 
contribution to the correlator (and, as a consequence, to susceptibilities) as the 
critical point is approached. 

In a non-interacting gas in thermal equilibrium the correlator (5-2) vanishes 
unless p = k and a = /?.*) The hadrons, however, are interacting. One can ask a 
question: what is the effect of the interaction on the correlator (5-2)? The answer can 
be found to leading order. 56 ) The contribution is proportional to the amplitude of 

We are not considering HBT correlations, which are a finite size effect. 
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the forward scattering A pk ^ pk of the particles with momenta p and k. This is easy to 
understand using the following argument. The amplitude of the forward scattering 
shifts the energy of the 2-particle state relative to the sum of single particle energies. 
The statistical weight of the two particle state is therefore changed relative to the 
product of the single-particle weights. The difference is the two-particle correlator: 



(n p n k ) - (n p )(n k ) = f p f k (e 0El - 1) » /„/*,(-/?#/) ~ f p f k (3A pk ^ pk . (5-4) 




where f p is the equilibrium distribution function and Ej is the interaction energy. 
The exact formula, obtained using diagrammatic analysis, 56 ) contains additional 
factors (1 + / p )(l + /fe), which can be understood as Bose enhancement (stimulated 
emission) factors (or, in the case of fermions, (1 — / p )(l — fk) ~ Pauli blocking). 

Near the critical point the most singular contribution comes from the exchange 
of the sigma field quanta in the t channel. *) Since, by kinematics, the quanta carry 
zero momentum, the singular contribution is proportional to l/m 2 , which equals £ 2 
- the square of the sigma field correlation length. 

The absolute strength of the singu- 
larity depends on the coupling of the 
critical mode sigma to the correspond- 
ing hadron in Fig. 10, which is difficult 
to estimate reliably. Order of magni- 
tude estimates have been made for cou- 
pling to pions 42 ) and to protons 57 ^ . 

As an example of the singular con- 
tribution in Fig. 10 consider baryon 
number susceptibility. Let Q in equa- 
tion (5-3) be the baryon number. Then 
one can see that the 1 / to 2 , or £ 2 , singu- 
larity from Fig. 10 for scattering two baryons results in the divergence of the baryon 
number susceptibility (2-2) (77 = at this order). If only charged baryons are de- 
tected, the total baryon number cannot be measured event by event, but the number 
of protons is measurable. Since, according to Fig. 10, the proton number fluctua- 
tions should also be singular at the critical point, measurement of such fluctuations 
may provide a signal of the critical point. 57 ) 

In principle, knowing the correlator (5-2) one could make quantitative predic- 
tions for fluctuation measures used in experiment. In practice, calculating the cor- 
relator is a very difficult task (what interactions should be included and what is 
their strength?). Non-equilibrium effects make this task even more difficult. Near 
the critical point these complications become less relevant since, as long as we limit 
ourselves to the singular effects, we only need to consider contributions such as in 
Fig. 10. 



Fig. 10. Diagrammatic representation of the 
singular contribution to the correlator 
(An p An k )- 



*-* Strictly speaking, what we call here, for simplicity, "sigma" is a mixture (a linear combination) 
of chiral condensate, baryon density and energy density fluctuations. 
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5.2. Fluctuations, correlations, and acceptance 

Cumulative measures of fluctuations are often used to represent experimental 
results. These measures suffer an important drawback - they depend on the size and 
shape of the acceptance window of the detector. This makes comparison of different 
experiments, as well as an experiment to a theory, difficult. However, knowing 
certain properties of the correlator (5-2), it is possible to correct for acceptance in 
such comparisons. 

As an illustration consider event-by-event fluctuations of the mean transverse 
momentum px per particle. Most commonly used fluctuation measures are based on 
the width of the distribution of the event mean pj>, 0" 2 bc . Similar to (5-3), <7 2 be can 
be expressed through the correlator (5-2) : 42 ) 

*ebe = E A VT&k T (An p An k ), (5-5) 

where Apx = Pt ~ Pf an d (N) is the average multiplicity of accepted particles. In 
the thermodynamic limit (N) — > oo the fluctuation o" 2 be vanishes as 1/(N). Thus in 
this limit, the quantity (N)a^ he does not depend on the size of the system (N) and 
is therefore a natural subject of theoretical predictions. 

To make closer comparison to experiment, it is useful to exclude the diagonal 
terms p = k from the sum in (5-5), since they give the trivial statistical contribution 
(N)~ 1 af nc , where <7i nc is the r.m.s. width of the inclusive distribution of pr- The 
remaining off-diagonal terms in (5-5) give the nontrivial "dynamical fluctuation", 
experimentally obtained after the subtraction: 

<4n = "ebe " W-V? nc . (5-6) 

In an experiment, the sum in (5-5) is limited to p and k which fall within detector 
acceptance. Assume, for clarity, that the acceptance is limited in rapidity, i.e., y p 
and yk belong to an interval [y m in> Umax]- The cumulative measure cr 2 be , or cr^ , 
will then depend on y acc = y max — ymin- This dependence simplifies in two regimes 
of y acc . The boundary between the two regimes is determined by the characteristic 
range y COTT of the rapidity correlator of the fluctuations: 



dyn 



y&cc 



0(y acc ), for y acc < y C or r ; 

(5-7) 

( N ) a dvn ' for ^cc > 2/corr- 



In other words, cumulative measure (A r )cr 2 yn grows linearly with y acc for small accep- 
tance windows and saturates at its thermodynamic limit value when the acceptance 
range exceeds the correlation range. In most current experiments, the width of the 
rapidity window y acc is smaller or at most comparable to the typical range of the 
rapidity correlator y corr ~ 1. This means that in a typical experiment, for a cumula- 
tive measure, normalized to be finite in the thermodynamic limit, the experimentally 
observed magnitude is roughly proportional to the acceptance window size. 52 )' 58 )' 59 ) 
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5.3. Experimental results and concluding remarks 

As an example of the QCD phase diagram scan, the plot in Fig. 11 shows the 
results of the measurements of the pr fluctuations using a cumulative measure 2Jp T 
described in Ref. 60). No clear non-monotonous signal, which one would expect if 
the experiments probed the vicinity of the critical point, was found. 

It is also interesting to compare the magnitude of the observed fluctuations to 
the singular contribution expected near the critical point. 42 ) After correcting for 
acceptance using the method outlined in the previous section one finds: 

where G is the magnitude of the owi: coupling in the diagram Fig. 10 and £ = l/m a . 

It is important to note, that ob- 
servation of a large magnitude of fluc- 
tuations would not by itself constitute 
the signal of the critical point. There 
are many possible contributions to the 
fluctuations, which are difficult to esti- 
mate. The distinct signature of the crit- 
ical point is the non-monotonous behav- 
ior of fluctuation observables. 

Experiments at other energies, at 
CERN SPS, RHIC, and future GSI fa- 
cility, will be able to provide a com- 
plete scan of the reachable domain on 
the QCD phase diagram and either dis- 
cover or rule out the presence of the crit- 
ical point in this domain. 

This review focused mainly on the 
signatures of the QCD critical point 
based on the event-by-event fluctua- 
tions. Further study of the properties of 
the critical point may reveal other, perhaps, even more sensitive and experimentally 
cleaner signatures. For example, real-time correlation functions and non-equilibrium 
dynamics near the critical point deserve further investigation. 62 ) 

Finally, the lack of a controllable and reliable theoretical method to calculate 
coordinates of the critical point impairs our ability to perform a more focused search. 
It is hard to overemphasize the importance of such a theoretical method. 
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